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^1 \ We study the computational complexity of approximately counting the number of independent 

sets of a graph with maximum degree A. More generally, for an input graph G — (V,E) and an 
activity A > 0, we are interested in the quantity Zq(X) defined as the sum over independent sets 
/ weighted as w(I) = A' 7 . In statistical physics, Zq(\) is the partition function for the hard-core 
model, which is an idealized model of a gas where the particles have non-negligible size. Recently, 
\mJ • an interesting phase transition was shown to occur for the complexity of approximating the 

[ partition function. Weitz showed an FPAS for the partition function for any graph of maximum 

degree A when A is constant and A < A C (T A ) := (A - 1) A_1 /(A - 2) A . The quantity A C (T A ) 
O | is the critical point for the so-called uniqueness threshold on the infinite, regular tree of degree 

A. On the other side, Sly proved that there does not exist efficient (randomized) approximation 
algorithms for A C (T A ) < A < A C (T A ) + e(A), unless NP=RP, for some function e(A) > 0. We 
remove the upper bound in the assumptions of Sly's result for A ^ 4, 5, that is, we show that there 
does not exist efficient randomized approximation algorithms for all A > A C (T A ) for A = 3 and 
A > 6. Sly's inapproximability result uses a clever reduction, combined with a second-moment 
ly-^ | analysis of Mosscl, Weitz and Wormald which prove torpid mixing of the Glauber dynamics for 

■ sampling from the associated Gibbs distribution on almost every regular graph of degree A for 

the same range of A as in Sly's result. We extend Sly's result by improving upon the technical 
work of Mossel et al., via a more detailed analysis of independent sets in random regular graphs. 
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1 Introduction 

For a graph G = (V, E) and activity A > 0, the hard-core model is defined on the set 1(G) of 
independent sets of G where set I 6 1(G) has weight w(I) := A^L The so-called partition function 
for the model is defined as: 

Z G (X):= £ w(I)= £ AH. 
IeX(G) I€I(G) 

The Gibbs distribution [i is over the set 1(G) where fi(I) = w(I)/Zc(X). The case A = 1 is 
especially interesting from a combinatorial perspective, since the partition function is the number of 
independent sets in G and the Gibbs distribution is uniformly distributed over the set of independent 
sets. 

The hard-core model has received considerable attention in several fields. In statistical physics, 
it is studied as an idealized model of a gas where the gas particles have non-negligible size so 
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neighboring sites cannot simultaneously be occupied [3, 1]. The activity A corresponds to the 
fugacity of the gas. The model also arose in operations research in the study of communication 
networks [7]. 

We study the computational complexity of approximating the partition function. Valiant [15] 
proved that exactly computing the number of independent sets of an input graph G = (V, E) is 
^P-complete. Greenhill [4] proved that even when the input is restricted to graphs with maximum 
degree 3, it is still #P-complete. Hence, our focus is on approximating the partition function. 

Weitz [16] gave an FPAS (fully polynomial-time approximation scheme) for the partition function 
of graphs with maximum degree A when A is constant and A < A c (Ta) := (A — 1) A_1 /(A — 2) A . The 
activity A C (T A ) is the critical activity for the threshold of uniqueness/non-uniqueness of the infinite- 
volume Gibbs measures on the infinite A-regular tree [7]. Recently, Sly [11] proved that, unless 
NP = RP, for every A > 3, there exists a function e(A) > such that for graphs with maximum 
degree A there does not exist an FPRAS (fully-polynomial time randomized approximation scheme) 
for the partition function at activity A satisfying: 



It was conjectured in Sly [11] and Mossel et al. [10] that the inapproximability result holds for all 
A > A c (Ta). We almost resolve this conjecture, that is we prove the conjecture for all A with the 
exception of A G {4,5}. 

Theorem 1. Unless NP=RP, there does not exist an FPRAS for the partition function of the 
hard-core model for graphs of maximum degree at most A at activity A when: 

• A = 3 and A > A C (T 3 ) = 4; or 

• A > 6 and A > A c (Ta); or 

• A = 4 and A G (A C (T 4 ) = 1.6875,2.01538] U (4,+oo); or 

• A = 5 and AG (A C (T 5 ) = 256/243, 1.45641] U (1.6875, 2.01538] U (4, +oo). 

Sly's work utilizes earlier work of Mossel et al. [10] which studied the Glauber dynamics. The 
Glauber dynamics is a simple Markov chain (Xt) that is used to sample from the Gibbs distribution 
(and hence to approximate the partition function via now standard techniques, see [6, 12]). For 
an input graph G = (V,E) and activity A > 0, the state space of the chain is 1(G). From a state 
X t G 1(G), the transitions X t — > X t +i are defined by the following stochastic process: 

• Choose a vertex v uniformly at random from V. 



• If X' G 1(G), then set AQ_|_i = X' , otherwise set Xt+i = Xt. 

It is straightforward to verify that the Glauber dynamics is ergodic, and the unique stationary 
distribution is the Gibbs distribution. The mixing time T m j x is the minimum number of steps T 
from the worst initial state Xq, so that the distribution of Xt is within (total) variation distance 
< 1/4 of the stationary distribution. The chain is said to be rapidly mixing if the mixing time is 
polynomial in n = \V\, and it is said to be torpidly mixing if the mixing time is exponential in n 



A C (T A ) < A< A c (T A ) + e(A). 



• Let 




X t U {v} with probability A/(l + A) 
X t \ {v} with probability 1/(1 + A). 
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(for the purposes of this paper, that means 2~mi x — Gx.~p(^l(n))) . We refer the reader to Levin et al. 
[8] for a more thorough introduction to the Glauber dynamics. 

Mossel et al. [10] proved that the Glauber dynamics is torpidly mixing, for all A > 3, for graphs 
with maximum degree A when A satisfies (*). This result improved upon earlier work of Dyer et al. 
[2] which held for larger A, but not down to the critical activity A C (T A ). The torpid mixing result of 
Mossel et al. [10] follows immediately (via a conductance argument) from their main result that for 
a random A-regular bipartite graph, for A satisfying (★), an independent set drawn from the Gibbs 
distribution is "unbalanced" with high probability. 

The proof of Mossel et al. [10] is a technically involved second moment calculation that Sly [11] 
calls a "technical tour de force". Our main contribution is to improve upon Mossel et al.'s result, 
most notably, extending it to all A > A c (Ta) for A = 3. Our improved analysis comes from using a 
slightly different parameterization of the second moment of the partition function, which brings in 
symmetry, and allows for simpler proofs. 

To formally state our results for independent sets of random regular graphs, we need to partition 
the set of independent sets as follows. For a bipartite graph G = (V\ U V2, E) where |Vi| = | V2 1 = n, 
for S > 0, for i e {1,2}, let 

if = {I £ 1(G) : \InVi\ > \lnV 3 -i\+6n} 
denote the independent sets that are unbalanced and "biased" towards V{. Let 

i s B = {i e i(G) : \i n Vi\ < \i n y 3 -i| + Sn} 

denote the set of nearly balanced independent sets. 

Let G(n, A) denote the probability distribution over bipartite graphs with n + n vertices formed 
by taking the union of A random perfect matchings. Strictly speaking, this distribution is over 
multi-graphs. However, for independent sets the multiplicity of an edge does not matter so we can 
view Q(n, A) as a distribution over simple graphs with maximum degree A. Moreover, since our 
results hold asymptotically almost surely (a.a.s.) over Q(n, A), as noted in [10, Section 2.1], by 
standard arguments (see the note after the proof of [9, Theorem 4]), our results also hold a.a.s. for 
the uniform distribution over bipartite A-regular graphs. 

Theorem 2. For A = 3 and any A > A C (T A ), there exist a > 1 and 5 > such that, asymptotically 
almost surely, for a graph G chosen from Q(n, A), the Gibbs distribution ji satisfies: 

M4)<a-"min{MXf),^2)}- 

Therefore, the Glauber dynamics is torpidly mixing. 

This proves Conjecture 2.4 of [10] for the case d = 3. For A > 4 we can extend Mossel et al.'s 
results to a larger range of A than previously known. 

Theorem 3. For all A > 4 there exists e(A) > 0, for any A where A C (T A ) < A < A C (T A ) + e(A), 
there exist a > 1 and 5 > such that, asymptotically almost surely, for a graph G chosen from 
Q(n,A), the Gibbs distribution fi satisfies: 

MX^^a-^minMXf),^)}. 

Therefore, the Glauber dynamics is torpidly mixing. The function e(A) satisfies: 
• If A > 6, e(A) > A C (T A ) - A C (T A+1 ). 
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• For A = 5, e(5) > .402912. 

• For A = 4, e(4) > .327887. 

In the following section we look at the first and second moments of the partition function. We 
then state two technical lemmas (Lemma 10 and Lemma 11) from which Theorems 1, 2 and 3 easily 
follow using work of Sly [11] and Mossel et al. [10]. In Section 4 we prove the technical lemmas. Some 
of our proofs use Mathematica to prove inequalities involving rational functions, this is discussed in 
Section 2.5. 

2 Overview 

Proceeding as in [10] and [11], roughly speaking, to prove Theorems 1, 2 and 3 we need to prove that 
there exist graphs G whose partition function is close to the expected value (where the expectation 
is over a random G ~ Q(n, A)). At the heart of the argument lies a careful analysis of the first two 
moments of the partition function. The aim of this Section is to give a brief technical overview of 
the analysis. For more details, the reader is referred to [10] and [11]. 

In Section 2.1, we define the quantities which will be prevalent throughout the text. In Sec- 
tions 2.2 and 2.3, we revisit the first and second moments and state our main technical Lemma. 
We then prove Theorems 1, 2 and 3 in Section 2.4. In Section 2.5, we clarify our use of computer 
assistance for the proofs of some technical inequalities. 

2.1 Phase Transition Revisited 

Recall, for the infinite A-regular tree Ta, Kelly [7] showed that there is a phase transition at 
A C (T A ) = (A - 1) A ~V( A - 2 ) A - Formally, this phase transition can be defined in the following 
manner. Let Ti denote the complete tree of degree A and containing i levels. Let pi denote the 
marginal probability that the root is occupied in the Gibbs distribution on Xg. Note, for even £ 
the root is in the maximum independent set, whereas for odd t the root is not in the maximum 
independent set. Our interest is in comparing the marginal probability for the root in even versus 
odd sized trees. Hence, let 

p + = lim p2e and p~ = lim P2t+i- 

l-^oc I— >oo 

One can prove these limits exist by analyzing appropriate recurrences. The phase transition on the 
tree Ta captures whether p + equals (or not) p~ . For all A < A c (Ta), we have p + = p~, and let 
p* := p + = p~ . On the other side, for all A > A c (Ta), we have p + > p~. Mossel et al. [10] exhibited 
the critical role these quantities p + and p~ play in the analysis of the Gibbs distribution on random 
graphs G(n, A). 

2.2 First Moment of the Partition Function 

Let G ~ G(n, A). For a, j3 > 0, let 

= {I eZ G \\IC\Vi\ = an, \lnV 2 \ =/5n}, 

that is, a is the fraction of the vertices in V\ that are in the independent set and f3 is the fraction of 
the vertices in V2 for an independent set of a bipartite graph G. We consider only values of (a,P) 
in the region 

K = {(a,P)\a,p > and q + /3 < 1}, 
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since it is straightforward to see that for a graph G ~ Q(n,A), it deterministically holds that 
Xq 13 = whenever (a, (3) ^ TZ. For (a, f3) G TZ, define also 

i&x^ 

i.e., Zq is the total weight of independent sets in Xq 13 . The first moment of Zq is 
M Z°/\ = A<-«»(;„) (;„) (9p) A - e X p(„* l(a ,«), 



where 



( ) 



and H(x) = — rrlnx — (1 — :r)ln(l — x) is the entropy function. The asymptotic order of ~Kg[Z G '^\ 
follows easily by Stirling's approximation. 

The first moment was analyzed in the work of Dyer et al. [2]. We use the following lemma 
from Mossel et al. [10] that relates the properties of the first moment to p* ,p + and p~ . The most 
important aspect of this lemma is that in the non-uniqueness region ~¥ j q\Z C q^\ is maximized when 
a ^ f3 ([2]) and more specifically for (a, (3) = (p + ,p~) ([10]). 

Lemma 4 (Lemma 3.2 and Proposition 4.1 of Mossel et al. [10]). The following holds: 

1. the stationary point (a, [3) of <3?i over TZ is the solution to (3 = 4>{a) and a = 4>{j3), where 

and the solutions are exactly (p + ,p~), (p~ ,p + ), and (p* ,p*) when A > A c (Ta), and the unique 
solution is (p*,p*) when A < A c (Ta); 

2. when A < A c (Ta), (p* ,p*) is the unique maximum of $i over TZ, and when A > A c (Ta) ; 
(p + ,p~) and {p~,p + ) are the maxima o/$i overlZ, and (p*,p*) is not a local maximum; 

3. all local maxima of <3>i satisfy a + f3 + A(A — 2)a(3 < 1; 

4- p + ,p~ ,p* satisfy p~ < p* < p + and when A — > A c (Ta) from above, we have p*,p~ ,p + — > 1/A. 

For every A > 3, define the region 

TZ A = {(a, P) | a, (3 > and a + j3 + A(A - 2)a/3 < 1}. 

Part 3 of Lemma 4 establishes that the local maxima of $i lie in TZa- Note that for A > 3, we have 
TZa C TZ. Hence, the local maxima for all A > 3 lie in the interior of TZ. 

2.3 Second Moment of the Partition Function 

The second moment of Z^ satisfies ([10]) 

E g [(Z^) 2 } « exp(n-max$ 2 (a,/3,7,5,£)), 

7,<5,e 
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where 



$ 2 (a, P, 7, 5, e) = 2(a + p) ln(A) + H(a) + #1(7, a) + H x {a - 7, 1 - a) + #(/3) + /3) 
+#x(/3 - 5, 1 - /3) + A (#1(7, 1 - 2/3 + 5) - iT( 7 ) + fl^e, 1 - 2/3 + 5 - 7) 
+#i(a - 7 - e,/3 - <5) - Hi (a - 7, 1 - 7) + H^a - 7, 1 - (5 - 7 - e) 
-F 1 (a- 7 ,1-Q ! )), (2) 

and = —x ln(x) — (1 — x) ln(l — x), H\(x, y) = — x(ln(x) — ln(y)) + (x — y)(ln(y — x) — ln(y)). 

To make $2 well defined, the variables have to satisfy (a, f3) £ 1Z and 

7,<5,e>0, a-7-e>0, /3 - S > 0, l-2/3 + <5-7-e>0, 
l-o-/3-e>0, /3-<5 + e + 7-a>0. 

Lemma 4 tells us that in the non-uniqueness region (which is the region of interest in this paper) 
the first moment is maximized when (a,/3) is (p + ,p~) (or symmetrically, (p~,p + ))- To show that 
these unbiased configurations dominate the Gibbs distribution with high probability (as desired for 
Theorems 2 and 3 we will apply the second moment method, as used in [10]. 

To that end, we need to analyze the second moment for (a,/3) = (p + ,p~), and show that 
K g [{Z^) 2 ] = 0({Eg[Z%> 3 }) 2 ). To do that we need to show that for for (a, /3) = (p + ,p~),K g [(Z^) 2 ] 
is roughly determined by uncorrelated pairs of configurations. This crux of this is to show that $2 
is maximized when 7 = a 2 and 5 = /3 2 , which is detailed in the upcoming condition which was 
proposed in [11]. 

Condition 1 (Condition 1.2 1 of Sly [11]). There exists a constant x > such that when \p + — a\ < \ 
and \p~ — /3| < x then g a ^{^,5,e) := $ 2 (a, /3, 7, 5, e) achieves its unique maximum in the region (3) 
at the point 

( 7 VV*) = (a 2 ,/3 2 ,a(l-a-/3)). 

As mentioned earlier, Condition 1 implies that for (a, (3) = (p + ,p~), Eg[(Z^) 2 ] = 0((Eg[Z^^]) 2 ) . 
While the implicit constant in the latter equality is bigger than one, an application of the small graph 
conditioning method [5, 17] shows that for (a, (3) = (p + ,p~), Zq is concentrated around its ex- 
pected value (up to a multiplicative arbitrarily small polynomial factor). This yields a lower bound 
on the partition function Zq for (a,/3) = (p + ,p~). On the other hand, it is straightforward to 
an upper bound on the partition function for balanced configurations a = f3. Consequently, one 
obtains Theorems 2 and 3 as detailed below in Section 2.4, which implies that the Gibbs distribution 
is unbalanced with high probability. Sly [11] uses random regular bipartite graphs as a gadget in 
his reduction and utilizes this bimodality property of the Gibbs distribution. 

Before stating our new results on when Condition 1 holds, it is useful to remind the reader the 
previously known values of A for which Condition 1 holds: 

• A > 3, and A c (Ta) < A < A c (Ta)+£(A) for some (small) e(A) > 0, ([10, Lemma 6.10, Lemma 
5.1]); 

• A = 6 and A = 1, ([11, Section 5]). 

Let Ai/ 2 (Ta) be the smallest value of A such that 0(0(1/2)) = 1/2 (0 is the function defined in 
Lemma 4). Equivalently A^C^a) is the minimum solution of 

/ 1 . \ l-l/A / , . . \ l/A 

(l + (l/A)VA) (l-(l/A)VA) =1. (4) 

1 The numbering in this paper for results from Sly's work [11] refer to the arXiv version of his paper. 
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The following Lemma is the technical core of this work. 
Lemma 5. Condition 1 holds for 

1. A = 3 and A > A C (T A ), and 

2. A > 3 and A G (A C (T A ), A 1/2 (T A )], 

Lemma 5 is proved in Section 4. As a corollary of Lemma 5 we get that Condition 1 holds for 
the range of A specified in Theorems 2 and 3. 

Corollary 6. Condition 1 holds for: 

1. For A = 3 and A > A C (T 3 ). 

2. For A > 6 and A C (T A ) < A < Ai /2 (T A ) and A 1/2 (T A ) > A C (T A _ X ). 

3. For A = 5 and A C (T 5 ) < A < A C (T 5 ) + .402912. 

4. For A = 4, A C (T 4 ) < A < A C (T 4 ) + .327887. 

Proof. Part 1 is identical to the first bullet in Lemma 5. 

Part 2 follows from the second bullet of Lemma 5 and the fact that for A > 6, it holds A 1 / 2 (T A ) > 
A c (T A _i). To see this, by (4), we have that 1 — (1/A 1 / 2 (T A )) 1 / A > 0, which implies that Ai/ 2 (T A ) > 
1. For A > 6, we have A C (T A ) < 1. Hence, for A > 7, we have A!/ 2 (T A ) > A c (T A _i). For A = 6, 
the claim follows from the fact that A C (T 5 ) = 256/243 < A 1/2 (T 6 ) » 1.23105. 

For A = 5, note that 

Ai/ 2 (Tb) - A C (T 5 ) > 1.45641 - 256/243 > .402912, 

which proves Part 3. 
For A = 4, note that 

A1/2 (T 4 ) - A C (T 4 ) > 2.015387 - 27/16 = .327887, 
which proves Part 4. □ 

Theorems 2 and 3 now follow from Corollary 6 as outlined earlier, and detailed in the following 
subsection. 

2.4 Proofs of Main Theorems 

We now proceed to prove Theorems 2 and 3. 

Proofs of Theorems 2 and 3. The proof is essentially the same as the proof of [10, Theorem 2.2] 
with minor modifications. We include the proof for the sake of completeness. 
Choose 5 > such that for 

X := min $i(x,y), Y := max $i{x,y) 

\x—p+\<5, \y— p~ \<o \x— y\<o 

it holds that r := X — Y > 0. To see that this is possible, note that $i is continuous and hence it 
is uniformly continuous at any closed and bounded region. Since $i exhibits a global maximum at 
(p + ,p~), the existence of 5 follows. 
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By Markov's inequality, we obtain that a.a.s. 



Zg Zg 

To bound mm{n(lf), we need the following result from [10]. While their result is only 

stated for (a, (3) close to (1/A,1/A), it can readily be verified (as Sly also observed, e.g., see the 
discussion before Lemma 3.4 in [11]) that their proof holds in a neighbourhood of {p + ,p~), whenever 
Condition 1 holds. 

Lemma 7 (Theorem 3.4 of Mossel et al. [10]). 2 Let A > 3. Suppose that Condition 1 holds, then 
for all sufficiently large n there exist (a n ,f3 n ) where 

a n = p + + o(l), j3 n = p~ + o(l), and na n and nj3 n are integers, (6) 

and it holds a.a.s. that 

For all n large enough, there exist (a n , f3 n ) where \a n —p + \ < 1/ra, |/3„ — p~\ < 1/n and na n and 
nf3 n are integers, and therefore (6) holds. 

From Lemma 7 and Corollary 6 it follows that a.a.s. Zq^ > exp(n(X — j)) and consequently 

, exp(n(X - f )) 



MP?) > 



Since both <]?i and $2 are symmetric with respect to a, f3, a similar statement to Lemma 7 holds 
with the roles of p + ,p~ interchanged, so that 



. s s exp(n(X - |)) 
mm{^(X 1 ),/i(X 2 )} > — 2_. (7) 

Combining (5) and (7), we obtain 

M4) < exp(-^) min{/i(Xf ), M (X|)}. (8) 

This completes the proof. 

The torpid mixing of the Glauber dynamics claimed in Theorems 2 and 3 follows from (8) by 
Claim 2.3 in [2], which is a standard conductance argument. □ 

For Theorem 1, we will use Lemma 5 combined with the work of Sly [11], but we need one ad- 
ditional ingredient. The following combinatorial result will be used to extend the inapproximability 
result for A = 3 to a range of A for A > 6. 

Lemma 8. Let G be a graph of maximum degree A and let k > 1 be an integer. Consider the graph 
H obtained from G be replacing each vertex by k copies of that vertex and each edge by the complete 
bipartite graph between the corresponding copies. Then, 



Z G ((1 + X) k -1) = Z H (X). 



2 The stated version of the theorem differs slightly from the version in [10]. In particular, in [10] (a,/3) is fixed, 
in the sense that it's independent of n. Here, (a.,/3) depends on n. In the application of this theorem in the proof 
of Theorem 2.2 in [10], it is unclear how they deduce the existence of a fixed (a, /3) and this is why we modified the 
statement of the theorem. Their proof of Theorem 3.4 in [10] still goes through for this slightly modified version. 
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Proof. Consider the map / : 1h — > %G that maps an independent set I of H to an independent set 
J of G such that v G J if and only if at least one of the k copies of v in H are in /. 

For an independent set J of G the total contribution of sets in /~ 1 (J) to Zh{\) is ((1 + A) fc — 1)' J ', 
since for each v G J we can choose any non-empty subset of its k copies in H. □ 

We are now ready to prove Theorem 1. 

Proof of Theorem 1. Sly's reduction [11, Section 2] establishes that for any A > A c (Ta) such that 
Condition 1 holds and also the following two inequalities hold: 

3 

(A - l)p + p~ < (1 -p + )(l -p~) and p + < -(l-p~), (9) 

5 

then there does not exist (assuming NP^RP) an FPRAS for the partition function of the hard-core 
model with activity A. For the two inequalities in (9), as Sly [11] points out, they are unnecessary. 
First off, the inequality (A — l)p + p~ < (1 — p + )(l — p~) is implied by Part 3 of Lemma 4. And the 
condition p + < |(1 — p~) is used to simplify the proof of Lemma 4.2 in [11]. For completeness, a 
slight modification of Sly's Lemma 4.2 without (9) in the hypothesis is proved in Section 3. 

Corollary 6 establishes that Condition 1 holds for a range of activities A for each A > 3. By 
the discussion above, we obtain that there does not exist (assuming NP^RP) an FPRAS for the 
partition function of the hard-core model with activity 

• AG (A c (T 3 ),oo] for A = 3. 

• AG (A c (T A ),A c (T A _i)] for A > 6. 

• AG (A C (T 5 ), A C (T 5 ) + .402912] for A = 5. 

• AG (A C (T 4 ), A C (T 4 ) + .327887] for A = 4. 

Lemma 8 (used with k = 2 and A = 3) implies that there does not exist (assuming NP^RP) 
an FPRAS for the partition function of the hard-core model with activity A > y/5 — 1 in graphs of 
maximum degree 6. 

The range of A for A and A — 1 where we can prove hardness (that is, (A c (Ta)j A c (Ta— i)]) 
overlap for A > 6. This is useful since the hardness for A — 1 automatically gives hardness for A. 
Thus for A > 6 we have the hardness result on the set 

(A C (T A ), Ac(Ta-i)] U • • • U (A c (T 6 ), A c (T 5 )] U (A c (T 5 ), A c (T 5 ) + 0.402912] U (V5- 1, oo) = (A C (T A ), oo). 

This concludes the proof of the theorem for A > 6. 

For A = 4, we have the hardness result on the set (A C (T4), 2.01538] U (4, +oo). 
For A = 5, we have the hardness result on the set 

(A C (T 5 ), 1.45641] U (1.6875, 2.01538] U (4, +oo). 

□ 

2.5 On the Use of Computational Assistance 

We use Mathematica to prove several inequalities involving rational functions in regions bounded by 
rational functions. Such inequalities are known to be decidable by Tarski's quantifier elimination [14], 
the particular version of Collins algebraic decomposition (CAD) used by Mathematica's Resolve 
command is described in [13]. The algorithms are guaranteed to return correct answers — they do 
not suffer from precision issues since they use interval arithmetic (a real number is represented using 
an interval whose endpoints are rational numbers). 
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3 Non-Reconstruction revisited 



In this section, we reprove Lemma 4.2 from Sly [11] without (9) in the hypothesis. This will allow 
us to focus on proving Lemma 5. 

Recall, Ta is the infinite A-regular tree, and p + ,p~ denote the marginal probabilities that the 
root is occupied for the limit of even and odd, respectively, sized trees. Let Ta denote the infinite 
(A — l)-ary tree rooted at p. (Thus, these two trees only differ at the root.) For fel, let Ti denote 
the tree with branching factor A — 1 and containing i levels. Let qi denote the marginal probability 
that the root is occupied in the Gibbs distribution on Tg. Analogous to p + and p~ , let 

q + = lim q 2 g and q~ = lim q<u+i- 

The densities q + ,q~ are related to p + ,p~ by: 

+ P + V 



1 — p 1 — p + 

There are two semi-translation invariant measures £i + and £i~ on Ta obtained by taking the 
weak limit of the hard-core measure of even-sized trees T21 and odd-sized trees T21+1, respectively. 
In these measures p + , jX~ on the infinite tree, q + and q~ , respectively, are the marginal probabilities 
that the root is occupied. These measures jx + and £i~ can also be generated by a broadcasting 
process, see [11, Section 4]. 

For v G Ta, denote by S Vi £ the vertices at level i in the subtree of Ta rooted at v. Let X p ^ t+ 
denote the marginal probability that the root p is occupied in an independent set X generated by 
the following process: we first sample an independent set X from the measure fi + , then we condition 
on the configuration X$ on S p /, and finally we sample an independent set X from the hard-core 
measure on Tg conditioned on X$ = Xg- Note that the configuration of the vertices in S p / is a 
random vector, so that X p ^^ + is a random variable. Define similarly X p ^_ for fi~ . 

We may extend this definition to an arbitrary vertex v E Ta at distance D from the root p, by 
setting X v ^ + = X p ^^ + if D is even or X v ^ + = X p ^_ if D is odd. Thus, X v ^^ + is the probability 
that, in an appropriate translation of p + , v is occupied conditioning on the configuration of S v n- 
Define similarly X v n~ for p~ . 

In Sly [11], it is proved that X Pi £ j+ (X p ^ _) is strongly concentrated around q + (q~ , respectively) 
under the condition that p + < |(1 — p~). This concentration is used by Sly for establishing the 
properties of the gadget he uses in his reduction. We need to reprove the concentration without the 
condition p + < |(1 — p~). 

Lemma 9 (see Lemma 4.2 in Sly [11]). When A > A c (Ta), there exist constants Ci(A, A), C2(A, A) > 
0, for all sufficiently large i, 



X P,i 



> exp(-Ci^) < exp(-exp(C20)- 

> exp(-Ci^) < exp(- exp(C20)- 



Proof. We follow the proof of [11]. For a vertex v G Ta, let N(v) = N 1 ^) denote its children, 
and for i 6 N, let N l (v) denote its descendants i levels below. For 1 < i < i and s £ {+,—}, let 
X N i^ t £_ itS denote the vector {X w £_ i s : w G N l (v)}. 

By standard tree recursions, it can be proved that for s G {+, — }, 

_ , An„ e]VW (i-^_i, s ) 

A «As — n ^N(v),i-l,s) 



1 + ^HweN^i 1 - X w,i-l,s) 
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In our framework, one may establish a contraction property for h, but it is slightly more straight- 
forward to look at depth two of the recursion. Note that 



l-X„ 



1 



i+An 



so that 



A,. 



w&N(v) 



1 



l-X, 



i+An 



weN{v) 1+A1L 6JV(M) (1-X V _ 2 , S ) 
A 



=: r(Af 7V 2( p )^_2, s )- 
By recursively applying (10), we obtain that 



(10) 



X 



g{ X N 2L (p)/-2L,s)- 



Note that r(X N 2 {p)/ _ 2 ^ s ) = g{X^{ P ),i-2,s)- 

For the rest of the proof, we focus on the case s = + and v being at even distance from the 
root p, the other cases being almost identical (since we are looking at depth two of the recursion) . 



Assume X, 



To prove the new version of the lemma, it suffices to obtain the following contraction property. We 
will prove for all £ sufficiently large, for all v £ lA, for all pairs X N 2^ ^_2,+> ^'n 2 {v) 1-2 + that differ 
at a single vertex z* , 



are two vectors which are equal except at one vertex z*GN (v). 



N 2 {v),1-2,+j 



< 



1 



+ 



1 



\x z * 



-2,+ 



(A-l) 2 ' 10 

To accomplish this, it suffices to prove that for £ large enough, it holds that: 



XL 



Or 



-(X N 2( v ) r 



< 



dX z *^2,+ 

Let w* £ N(v) be the parent of z*. Then (12) is equivalent to: 



1 1 
(A -l) 2 + 10' 



(11) 



(12) 



A 2 n 



z£N(w*)\{z*} 



'l-X 



zJ-2 



w£N(v)\{w*} 



i + An 



z£N(w) 



'l-X 



z,e-2,+. 



A+n 



i + xU 



z£N(w) 



,1 — X Zi £-2,+) 



< 



1 

+ "777 • 



(A-l) 2 10 



Note that as £ — > +oo, since v is at even distance from p, X V £ + converges almost surely to 
(see [11, equation (4.1)] and [11, proof of Lemma 4.2]). Hence, as £ — > +oo, we have that 



Or 



QX 7 , 



z*,e-2,+ 



A 2 (l- g +)*-*(l + A(l- 9 +) 



A-nA-2 



A + (l + A(l-g+) A - 1 ) A ~ 1 



=: 7. 



Pick Iq = £o(A, A) so that for all £ > £q we have that 

Or 



OX 



z*,i-2,+ 



(Xn 2 (v),i-2,+) ~ 7 



1 

< — . 
- 10 
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We will see next that 7 = q + q~~ and hence by Part 3 of Lemma 4, it follows that 7 < r%-Tp wn i cn 
proves (12). 

To see that 7 = q + q~ , note that A(l — q + ) A ~ 1 = 1 1 _ and A(l — g _ ) A ~ 1 = jz^t ( one can derive 
these equalities by p + = (fi(p~),p~ = <^>(p + ), see Part 1 of Lemma 4). It is now a matter of a few 
algebra substitutions to check that 7 = q + q~ ■ 

This proves (11). Recursively applying this relation (11) implies that for all i — L > Iq, for 
A > 3, 

/j \ LV2J 

\9( X NL(v),e-L,+) - 9(XnL(v),£-L,+ )\ - [2) \ x z*/-l,+ - X 'z*,e-L,+ \- ( 13 ) 

Having established (13), the rest of Sly's proof of Lemma 4.2 goes through. □ 

4 Analysis of the Second Moment 

In this section, we prove Lemma 5. We break its proof into two slightly easier components, namely 
Lemmata 10 and 11 stated below, depending on the value of p + . Notice that Lemma 10 determines 
explicitly a region of (a, ft) for which Condition 1 holds. 

Lemma 10. Let A > 3 and let (a,/3) S TZa, 01, ft > 0, and a, ft < 1/2. Then g a 3(7, 5, s) := 
^(Qo Pi 7) 8, e) achieves its unique maximum in the region (3) at the point (7*, <$*,£*). 

Lemma 11. Fix A = 3 and A > A c (Ta). Let p + and p~ be the corresponding probabilities. Assume 
that 1/2 < p + < 1. There exists a constant x > such that for \p + — a\ < x an d \p~ — P\ < X, 
g a 0(7, 5, e) := ^2(0, ft, 7, 5, e) achieves its unique maximum in the region (3) at the point (7*, 5*, £*). 

Proof of Lemma 5. We begin by proving the second part of the lemma (A > 3 case) and then we 
prove the first part (A = 3 case). 

Fix any A > 3 and A > A c (Ta). Let p + , p~ be the corresponding marginal probabilities 
that the root is occupied in the measures fj,~. By the third item of Lemma 4, we have that 
(a, ft) = {p + ,p~) is contained in the interior of Hence, the same holds for every (a, ft) in a 
small enough neighborhood of (p + ,p~)- Hence, provided that p + < 1/2, Lemma 10 verifies the 
condition for any such point. It is now easy to see that p + is increasing in A, and hence p + < 1/2 
iff A G (A c (Ta), A^C^a)]- This proves the second part of Lemma 5. 

The first part of Lemma 5 is proved analogously. When p + < 1/2, one uses Lemma 10 as above, 
while when p + > 1/2 the condition reduces to Lemma 11. □ 

Thus we may focus our attention on proving Lemmata 10 and 11. While the analysis at some 
point for Lemma 11 requires tighter arguments, the two proofs share many common preprocessing 
steps. The rest of this section is devoted to these common preprocessing steps and the proofs of 
Lemmata 10 and 11 are given in Section 5. 
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4.1 The Partial Derivatives 



The derivatives of $2 with respect to 7, 5, e can easily be computed and are also given in [10, Proof 
of Lemma 6.1]: 

(d$2\ (1- 2(3 + 5 - 7 - e) A (a -7- e) A (1 - 2a + 7 ) A ^ 

eXV \djJ (l-/3-7-e) A (/3-a + 7-<5 + £) A (a- 7 ) A -V 

}e*£\ ((3-a-5 + 1 + e) A (l-2P + 6) A -i 
6XP \05 J (l-2/3 + 5-7-e) A (/3-5) A - 2 5 ' 1 j 

(d$2\ (l-2/3 + ^-7-6) A (a-7-e) A (l-«-/3-g) A , . 

eXP \de) e A ((3-a-5 + 7 + e) A (l-(3- 7 -e) A ' 1 j 

4.2 Excluding the Boundary of the Region 

We argue that for (ct,(3) G 7£a 5 the maximum of g a ^ cannot occur on the boundary of the region 
defined by (3). 

Lemma 12. For every A > 3 and (a, (3) G 7t&, <?a,^(7, <5, e) '■= ^2(0, /3, 7, 5, e) attains its maximum 
in the interior of (3). 

Proof. We follow the proof of Lemma 6.1 of [10] (where the same result is proved in the special case 
a = (3 = ±). 

We will prove that g a #(7, <5, £) attains its maximum in the interior of (3) by showing that at 
least one derivative in (14)-(16) goes to infinity (+ or — according to the direction), as we approach 
one of the boundaries defined by (3) from the interior of (3). 

For A > 3, note that (a, 0) G 72-a implies a + (3 < 1. Without loss of generality, we assume 
that a > j3. Hence f3 < 1/2. We have (15) goes to +00 as S — > 0, and (16) goes to +00 as e — > 0. 
We also have (16) goes to —00 as 7 + e — > a, (15) goes to —00 as S — > j3, (16) goes to —00 as 
7 + e — 5 — > 1 — 2(3, (16) goes to —00 as e — > 1 — a — (3, and (16) goes to +00 as 7 + 6 — 5 — > a — (3. 

When a > 1/2, the condition 7 = is not a boundary, as 7 > is implied by the conditions 
5 > 0, 1 — a — f3 — e>0, and (3 — 5 + e + "f — a > 0. On the other hand, when a < 1/2, we have (14) 
goes to +00 as 7 — > 0. □ 

4.3 Eliminating one variable 

Fix A,q,/3,7, 5 and view $2 as a function of e. We maximize with respect to e. In this setting, it 
was proved [10, Lemma 6.3], that the only maximizer of the function $2 in the interior of (3) is 
obtained by solving = and is given by: 

e:=e{a,M,S) = \{l + a - P - 2 1 - ^D), (17) 



where 



Define 



D = (l + a-/3-2 7 ) 2 -4(a- 7 )(l - 2(3-^ + 5) 
= (a + /3-l) 2 + 4(a- 7 )(/3-5). 



f!:=f)(a,f3 n ,5) = ^(1 - a + - 26 - </D), (18) 
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and note that 



a _ 7 _e = /3 _ S -fj = -(-(1- a-fl + VB), 

(a — 7 — e)(l — a — ft — £ — fj) = irj. 

The new parameter fj (not used in [10]) is symmetric with e, i.e, the constraints and formulas we 
have are invariant under a symmetry that swaps a, 7, e with /3,S,fj. This will allow for simpler 
arguments (using the symmetry). 

From the previous discussion and equation (17), we may eliminate variable e of our consideration. 
Of course, this introduces some complexity due to the radical y/~D, but still this is manageable. Let 

/(7,<5) := g a>l 3(7,S,£) = <S> 2 (a,f3,-f,5,£). 

To prove that (7*, 5* ,£*) is the unique global maximum of g a> p in the interior of the region defined 
by (3), it suffices to prove that (7*, 8*) is the unique global maximum of / for (7, 5) in the interior 
of the following region, which contains the (7, <5)-projection of the region defined by (3): 

0<7<a, 0<<5</3, 0< 1-2/3 + 5-7, < 1 - 2a + 7 - 5. (20) 

Each inequality in (20) is implied by the inequalities in (3), the only non-trivial case being the last 
inequality which is the sum ofl — a — /3 — e > and (3 — 5 + e + j — a>0. 
The first derivatives of / with respect to 7, S are 

df 

7T-(7,<S) = AlnWii+lnWia, (21) 



df 

-^( 7 ,<5) = AlnW 2 i+ln^ 22 , (22) 
00 



where 



(q-7-g)£(l-2a + 7) e(l-2a + 7 ) (a- 7 ) 5 

vVll = —, tt\ = ~7Z n s") W 12 ~ 



rj(a — 7) 2 (1 — a — (3 — e)(a — 7) ' (1 — 20 + 7)7' 

(P-5-fj)f)(l-2p + 5) fj{\-2fi + 5) (/3-5) 2 

W 21 = T7T, FTo = T-, 7, ^T77^ FT' ^22 ~ 



e(/3-<5) 2 (1 -a-P-rj)(P-6y (1-2/3 + 5)5' 

Note that the rightmost equalities in the definition of W\\ and W%x follow from (18) and (19). 

For every A > 3, and (a,/3) £ TZa, we have that (7*, 5*) is a stationary point of / (this follows 
from the fact that for 7 = a 2 and S = /3 2 , the inequalities on the right-hand sides in Lemma 13 
become equalities, and from (21), (22) we have that the derivatives of / vanish). 

4.4 Restricting the Region 

To determine whether (21) and (22) are zero it will be useful to understand conditions that make 
Wix, W12, W21, W22 greater or equal to 1. The following lemma gives such conditions. The proof is 
given in Section 6.1. 

Lemma 13. For every (a, (3) € 1Z, and (7, <5) in the interior of (20), 

W U >1 (l-a) 2 <5 + /3 2 (2a-l- 7 ) > 0, 

W 12 > 1 7 < a 2 , 

W 2 i>l ^ (l-/3) 2 7 + a 2 (2/3-l-5) > 0, 

W 22 > 1 ^ 5 < /3 2 . 
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By considering the sign of (21) and (22) and Lemma 13, we have that the stationary points of 
/ (and hence of g a ,/s) can only be in 

< 7 < a 2 , < 5 < f3 2 , (1 - a) 2 5 + /3 2 (2a - 1 - 7) < 0, (1 - /3) 2 7 + a 2 (2/3 - 1 - 8) < 0, (23) 

or 

a 2 < 7 < a, P 2 < 5 < (3, (1 - a) 2 5 + (3 2 (2a - 1 - 7) > 0, (1 - /3) 2 ~f + a 2 (2/3 - 1 - 6) > 0. (24) 

Note that for 7 = a 2 , 5 = f3 2 one has W\\ = W12 = W21 = W22 = 1 (Lemma 13 holds with equalities 
as well instead of inequalities), so that (a 2 ,/3 2 ) is always a stationary point for 7(7, 5). 

4.5 The Hessian 

To prove that / has a unique maximum, we are going to argue that / is strictly concave in each of 
the regions defined by (23) and (24). It will thus be crucial to study the Hessian of /. 
Let H denote the Hessian of /, i.e., 



J§(7,*) B§fe(7,*)' 



H 



Our goal is to express H in a helpful explicit form. In this vein, it will be convenient to define 
the following quantities. 

R _ 1-0-/3 „ _ y/p „ _ 2(q-7-e) 

1X1 ~ l-u-P-e-r)i 1X2 ~ l-2a+7' ^3 — Q _ 7 , 

n _ \[D o _ 2(1-/3-7-e) „ _ y/p 

^4 - — , H 5 - — , Rq - jz^g+y, 

_ 2(a- 7 -e) R _ o _ 2(l-/3- 7 -e) 
«7 - ' -"8 - — , K 9 - — 5 . 

H can now be written in a relatively nice form with respect to the Ri. Namely, 

df 1 r i 

-H 7 ,5) = -= + i? 2 + i? 3 )A - i? 2 - i? 3 - Ra - R 5 , (25) 

o z j s/D L -I 

<9f 1 r l 

^(Tf,S) = -= (-Ri + R 6 + i? 7 )A - R 6 - R 7 - R 8 - R 9 , (26) 

dW5^ 5) = 9^ (7 ' 5) = W (27) 

Inspecting the i?j we obtain the following observation. 
Observation 14. Ri, . . . , Rg are positive when (a,/3) G 7£a (7><5) * n ^ e interior of (20). 

Observation 14 and equation (27) immediately yield: 
Observation 15. For every (a,/3) G T^Aj and (7,5) in the interior of (20), 

f :(7,*) = Jln(7» > 0. 



In Section 6.2 we prove the following technical inequality on the Ri. 
Lemma 16. For every (a,/3) £ 1Z, and (7,(5) in the interior of (20), 

Ri> R 2 + i?3, and R\ > R & + R-j. 
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Applying Lemma 16 and Observation 14, (25) and (26) give the following straightforward corol- 
lary. 

Corollary 17. For every (a, (3) 6 TZ, and (7,(5) in the interior of (20), 

J^XO, |(7^)<0. 

Corollary 17 implies that the sum of the eigenvalues of M is negative. Hence, H is negative 
definite iff the determinant of H is negative. We have the following expression for det(H). 



df df 



d f ( x\ °l r x\ 



1 

D 



Rj 



(A - If + R 2 + R 3 )(-Ri + Re + Rr) - 
+(A - 1) + R 2 + R 3 )(-Ri -Rs- Rg) + 

+ (—R\ + Re + Rt){—Rx — Ra — R5) — 2R\ 
{-Rx -R 8 - Rg)(-Ri -Ra- R 5 ) - R\ 



(28) 



5 Concluding the Proofs of Lemmata 10 and 11 

In this section, we give the proofs of Lemmata 10 and 11. We first recap what we have accomplished 
in Section 4 for every (a, (3) € TZa- 

1. The function <7 a 5(7, <5, e) attains its maximum in the interior of the region (3). See Section 4.2. 

2. To study the (local) maxima of g in the interior of the region (3), it suffices to study the maxima 
of the function 7(7, 5) = g a ^{^, 5, e) in the interior of the region (20). More explicitly, if 7(7, S) 
has a unique maximum in the interior of the region (20) at (7*, 5*) = (a 2 , /3 2 ), then g a p has a 
unique maximum in the interior of the region (3) at (7*, 5* , £*) = (a 2 , /3 2 , a(l — a — /3)) . The 
function / is differentiable in the interior of the region (20). See Section 4.3. 

3. The point (7*, 5*) = (a 2 , /3 2 ) is always stationary for /. Every stationary point of / lies in one 
of the two regions defined by (i) (20) and (23), (ii) (20) and (24). See Section 4.4. 

4. The function / is strictly concave iff det(H) > 0. See Section 4.5. 

By the above discussion, if for some (a, /?) £ TZa it holds that / is strictly concave in each of 
the regions (i) (20) and (23), (ii) (20) and (24), then g a ^ has a unique maximum at (7*, £*,£*) = 
(a 2 ,/3 2 ,a(l — a — /?)). Thus, it suffices to check that det(H) > in each of these regions. This is 
essentially the way we derive Lemmata 10 and 11. 

Hence, the main challenge is proving that det(H) > 0. This can be done slightly more easily in 
the region (20) and (23). Indeed, in Section 6.3 we prove the following lemma. 

Lemma 18. det(H) > for every A > 3, (a, {$) G TZa, (7, * n ^ e interior of (20) and (7,(5) 
in (23). 

Proving det(H) > in the intersection of the regions (20) and (23) is trickier and is essentially 
the reason we do not obtain our hardness result for A = 4, 5. At this point, it is convenient to split 
the analysis for each of the lemmas. 
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5.1 Proof of Lemma 10 

In the setup of Lemma 10, we have (a, f3) € an d a,/3< 1/2. We suppress the details of proving 
det(H) > as a lemma, whose proof we defer to Section 6.4. 

Lemma 19. det(H) > for every A > 3, (a, (3) G Ha> ot, /3 < 1/2, (7,5) in the interior of (20) 
and (7, 5) in (24)- 

Using Lemmata 18 and 19, the proof of Lemma 10 is immediate. 

Proof of Lemma 10. Lemma 18 and Lemma 19 imply that / has a unique maximum at (7*, 5*) for 
every A > 3, (a,/3) £ 1/2 > a, (3 and (7, 5) in the interior of (20). This follows from the fact 
that det(H) > implies that the Hessian of / is negative definite in the region of interest, i.e., where 
g a< p could possibly have stationary points, which in turn implies that / is strictly concave in the 
region and hence has a unique maximum. By the definition of /, it follows that g a a has a unique 
maximum at (7*, £*,£*). For a more thorough outline, see the beginning of Section 5. □ 



5.2 Proof of Lemma 11 

In the setup of Lemma 11, we have A = 3 and p + > 1/2. By (1), tedious but otherwise simple 
algebra gives that the solution of /3 = <j>(a) and a = <j)((3) with a 7^ /3 satisfies 

a 2 - 2a + a/3 + 1 - 2/3 + /3 2 = 0. 

It can also be checked that when 1 > p + > 1/2, it holds that < p _ < (3 — \/5)/4. Define 1Z' 3 to be 
the set of pairs (a,/3) such that a 2 -2a + a(3 + l-2(3 + f3 2 = 0, 1/2 < a < 1 and < f3 < (3-\/5)/4. 
Our goal is to show that det(H) > for every (a, (3) E 1Z' 3 , (7, 5) in the interior of (20) and (7,(5) 
in (24). 

We can rewrite det(H) using the formula in (28) as 

det(H) = ^\zRi{U 1 + U 2 ) + U 1 U2 



D 

Ui 
D 



3i?i(l + U 2 /U 1 ) + U 2 , (29) 



where 

JJ\ = Rg + Rg — 2Rq — 2Rj, U2 = R4 + R5 — 2R2 — 2i?3. 

The following lemma establishes technical inequalities on U\ , U2 , Ri , . . . , R9 , which are crucial 
in establishing the positivity of det(H). Its proof is given in Section 6.5. 

Lemma 20. We have U\ > 0, i? 4 > R 6 , R 5 > R 7 , R 5 > AR 3 , R 9 > AR 7 , and 3R 8 /2 + R 9 > 9R 2 , 
for every (a, (3) £ 1Z' 3 , (7, 5) in the interior of (20) and (7,(5) in (24)- 

With Lemma 20 at hand, we can now prove that det(H) is positive. 

Lemma 21. det(M) > for every (a, (3) G TZ' 3 , (7,(5) in the interior of (20) and (7,(5) in (24)- 

Proof. Observe that 

Ui + 3U 2 = (Rs + 2R 9 /3 - 6R 2 ) + (3i? 4 - 2R d ) + {R 9 /3 - 4i? 7 /3) 
+ (3i? 5 /2 - 2# 7 /2) + (3i? 5 /2 - 6R 3 ) 
>0, 
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where the last inequality follows by Lemma 21. Once again by Lemma 21, we have U\ > and 
hence U 2 /Ui > -1/3. Thus, (29) gives 

det(H) > ^(2i2i + C/ 2 ) = ^ [i? 4 + R 5 + 2(i?i - i? 2 - fls)] > 0, 

where the last inequality follows from Lemma 16 and Observation 14. □ 

Proof of Lemma 11. As in the proof of Lemma 10, Lemmata 18 and 21 yield that g a $ (7, S, e) 
achieves its unique maximum in the region (3) at the point (7*, £*,£*), for every (a, (3) S 1Z' 3 . We 
next show that g a ,p{l-> S, e) := ^2(0, [3, 7, 5, e) also achieves its unique maximum in the region (3) at 
the point (7* ,5* ,£*) in a small neighborhood of 1Z' 3 . 

First note that $2 is continuous. By Lemma 12, we have for sufficiently small x > 0, the 
maximum of g a ^ cannot be obtained on the boundary of the region (3). 

Note that the derivatives of $2 are continuous. It follows that for sufficiently small x > 0) 
all stationary points of g a ^ have to be close to the point (7*, £*,£■*). We can choose x such that 
det(H) > in the neighborhood of the point (7*, £*,£*), which implies that g a R has a unique 
stationary point and it is a maximum. □ 



6 Remaining Proofs of Technical Lemmas 
6.1 Proof of Lemma 13 

We define W 3 to be the numerator of W\\ minus the denominator of Wn, and W4 to be the numerator 
of W21 minus the denominator of W21', more precisely 

W 3 = (a-7-e)e(l-2a + 7)-r)(a-7) 2 , (30) 
W4, = W-5-f))f)(l-2P + 5)-e(f3-5) 2 . (31) 

Substituting the expression for i and simplifying we obtain 

W 3 := ((a + /3-l + \/^)/2)e(l-2a + 7)-r?(a-7) 2 , 
W 4 := ((a + P-l + VD)/2)f](l- 2/3 + 5) - e(0 - 5) 2 . 



By expanding W 3 and W4, we have 



where 



W 3 = 1^31+^32^, 

W 4 = Wu + W A2 ^D, 



W 31 = ((3/2)0 - (l/2)0 2 - 5 - (3/2)a0 + Sofa 

+0 + (5/2)0a 2 + (l/2)a 3 + (3/2)a - (3/2)a 2 
-(l/2)0 2 - 1/2 + 5a - (7/2)a0 - 5a 2 + 2 a, 

W 32 = a/3 - (1/2)07 " a " (V2)/3 + (l/2)a 2 + 1/2, 

W41 = ((3/2)a-(l/2)a 2 -7-(3/2)a/3 + 7 /3),5 

+a + (5/2)a0 2 + (l/2)0 3 + (3/2)0 - (3/2)0 2 
-(l/2)a 2 - 1/2 + 70 - (7/2)q0 - 7 2 + a 2 0, 

W42 = a0-(l/2)a5-0-(l/2)a + (l/2)0 2 + l/2. 
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Note that for every (a, 0) € TZ, W 32 > when < 7 < a, and W42 > when < 5 < 0. To see 
W 32 > note that a/3 - (1/2)07 -a- (1/2)0 + (l/2)a 2 + 1/2 > a0 - (l/2)0a -a- (1/2)0 + 
(l/2)a 2 + 1/2 = (1/2)(1 - a)(l - a - 0) > 0, since a + < 1; inequality W42 > is the same (after 
renaming the variables). 

Note that W31 is a linear function in 7 and we have 

^1 = (3/2)0 - (l/2)0 2 - 5 - (3/2)a0 + 5a, 

Of) 

which is positive for all (a, 0) £ 7£ and < 5 < 0. To see this note 

(3/2)0 - (l/2)0 2 -5- (3/2)a0 + <fa = (1 - a)((3/2)0 - 5) - (l/2)0 2 
> (1 - a)(l/2)0 - (l/2)0 2 = (1/2)0(1 - a - 0) > 0. 

Moreover, we have W31 is negative when 7 = a, (a, /3) £ TZ and < 5 < (after substituting 7 = a 
into W31 we obtain -(1 - a)(l - a - 0) 2 /2 < 0). Hence, W31 < for all (a, 0) G 7£, < 7 < a and 
< 5 < 0. 

By the same proof, we can show that W41 < for all (a, 0) £ 7£, < 7 < a and < 5 < (note 
that W31 and W41 are the same after renaming the variables). 
Let 

W 5 = (W31/W32) 2 - -D, 
W 6 = (W41/W42) 2 - D. 

Note that the signs of W3 and W5 are opposite, and the signs of W4 and Wq are opposite. 
After substituting W31, W32, W41, W42 and simplifications, we obtain 



W 5 
W 6 



4(0 - 5)(q - 7 ) 2 ((1 - a) 2 5 + 2 (2a - 1 - 7)) 

(2a0 - 07 - 2a - + a 2 + l) 2 
4(0 - 5) 2 (a - 7 )((1 - 0) 2 7 + a 2 (20 - 1 - 5)) 



Hence, 

and 



(0 2 - 5a + 1 - 20 - a + 2a0) 2 
W u > 1 W s <0 ^> (1 - a) 2 5 + 2 (2a - 1 - 7) > 0, 



W 2 i > 1 ^ W 6 <0 ^ (1 - 0) 2 7 + a 2 (20 - 1 - 5) > 0. 
Finally, we analyze the conditions for W12 > 1 and W22 > 1- It is straightforward to see that 

W12 > 1 (a - 7 ) 2 > (1 - 2a + 7)7 7 < a 2 , 

and 

W 22 > 1 <^=> (0 - ,5) 2 > (1 - 20 + 5)5 ^ 5 < 2 . 

6.2 Proof of Lemma 16 

We have 

--Rl + i?2 + i?3 = (0 " 5) r - - + - - -V + : I + 



a — 7 — i fj i J £ 1 — 2a + 7 a — 7 
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Multiplying by the denominators we let 



Pi = + R 2 + R 3 )£fj(l -2a + 7)(a - 7) (a - 7 - e) = P n + P12VA 



where 



P n = 1 - (l/2)*a + Hap + (9/2)* 7 - 5/3 7 + 21a 7 /3* - (37/2)a 2 7 /3* - 5a 7 /3* 2 + (5/2)a 7 /3 2 * 
+87 2 a5/3 - 5a - 3/3 - * + 10a 2 + 3/3 2 + 35(3 - 10a 3 + 5a 4 - /3 3 + 15a7/3 - (13/2)a/3* 
-(27/2)a 7 * - (5/2)07* - 4a 2 /3* + (27/2)a 2 7 * - 15a 2 7 /3 - 24a 7 /3 2 - 8 7 2 /3* 
+(17/2)a/3 2 * - (5/2)/? V - 3a/3* 2 + 3/? 7 * 2 + (15/2)a 3 /?* + 4a 2 ^<5 2 - (ll/2)a 2 /3 2 * 
-(9/2)a 3 7 * + 3a 2 7 * 2 + 5a 3 7 /3 + 17a 2 7 /3 2 + 4a 7 /? 3 - 3 7 2 a* 2 - 5 7 2 a/3 2 + 7 2 /3* 2 - 3/3 2 * 
+3* 2 a - 3* 2 7 - 3a 2 * 2 + (7/2)a 4 * - 8a 4 /? - 15a 3 /3 2 - 4a 2 /3 3 + 3 7 2 * 2 - 7 2 /? 3 - 33a 2 
+(15/2)a 2 * - 16a/3 2 + 7/? 2 7 - (19/2)a 3 * + 27a 3 (3 + 28a 2 /3 2 + 5 7 2 /3 2 + 4a/3 3 - 2/3 3 7 
+5(3 3 - (3/2)a/3 3 * + (l/2)/3 3 7 * - a 5 , 

and 

P 12 = _i_(i/2) ( 5a-9a/3-(5/2)*7 + 3 ) 57-2a7/3(5 + 4a + 2/3 + (5-6a 2 -/3 2 -2*/3 
+4a 3 - a 4 - 6a7/3 + 4a/3* + 5a7* - a 2 /?* - (5/2)aV + 3a 2 7/3 + 4aj(3 2 + 7 2 0* 
-(3/2)a,S 2 * + (l/2)/?V - «7* 2 + /3 2 * - * 2 a + * 2 7 + a 2 * 2 + 12a 2 /3 - 2a 2 * + 4a/3 2 
-2/3 2 7 + (3/2)a 3 * - 5a 3 /3 - 4a 2 /3 2 - 7 2 /3 2 . 

The following claim is proved using the Resolve function of the Mathematica system in Ap- 
pendix A. 2. 

Claim 22. Pn > and P i2 < for all (a, (3) G 1Z and (7,*) in the interior of (20). 



From Claim 22 we have that Pn — P\2\D > and hence showing Pn +P\2\D < is equivalent 
to showing P n — P 2 2 D < which in turn is equivalent to showing that P2 := (P11/P12) 2 — D is 
negative. We have 



where 



pi 



P 21 = -8*a + 20a 7 /3* - 10a 2 7 /3* + 2a 7 /3* 2 + 12a 7 /3 2 * + 25 + 25 2 - 850 + 4/3 3 + 34a/3* 

-IO/37* - 44a 2 /3* + 8a/3 2 * - 4/3 2 7 * - Ua05 2 - 2/3 7 * 2 + 18a 3 06 + 6a 2 05 2 - Wa 2 2 5 
+9a 2 7* 2 - 16a7^ 3 - 18a 7 * 2 - 20 2 5 - 15* 2 a + 9* 2 7 + 24a 2 * 2 + 2a 4 * + 16a 2 /3 3 
+4 7 2 ( S 3 + 12a 2 * - 8a 3 * - 16a,S 3 + 8/3 3 7 - A50 3 + 6a,S 3 * - 2/3 3 7 * + 85 2 /3 + 20 2 5 2 
+8* 3 a - 4* 3 a 2 - ll* 2 a 3 - 35 2 a0 2 + * 2 7 /3 2 - 4* 7 2 /3 2 - 4* 3 . 

The following claim is proved using the Resolve function of the Mathematica system in Ap- 
pendix A. 3. 

Claim 23. P21 > for all (a, (3) £ 1Z and (7,*) in the interior of (20). 



From Claim 23 we have that P2 < which implies Pn + P\2\D < and this completes the 
proof of Lemma 16. 
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6.3 Proof of Lemma 18 

We will use the following technical lemma in the proof of Lemma 18. 

Lemma 24. For every (a, (5) £ 1Z, and (7,(5) in the interior of (20), we have 

(-Pi + R 2 + R 3 )(-Ri + Re + Rt) -Rj<0, (32) 
(-Pi + P 2 + Ra)(-Ri -Us- R9) ~ R\ > 0, (33) 

and 

{-Ri + P 6 + R 7 )(-Ri -R A - R5) - P 2 > 0. (34) 
Proof of Lemma 24- We first prove inequality (32). We have that 

(— R\ + R2 + Rs){— Ri + Rq + R7) = (Pi — i?2 ~~ P3XP1 ~ R§ ~ R7) < P?, 

where the last inequality uses Lemma 16 and the positivity of the Ri (Observation 14). This 
establishes (32). 

We next prove (33), noting that inequality (34) follows by an analogous argument. The left-hand 
side of (33) multiplied by the denominators and simplified is 

P5 = {{-Ri+R2 + Rz){-Ri-Rs-Rv)-Rl)o~{a-i)sf l {l-2a + 1 ){u- 1 -e)/^ 

= p 51 + p 52 Vd, 

where 

P51 = -l-4/37 2 5a + 35a-17a/3 - 357 + 5/?7-14a/?75 + lla 2 /3 7 5 + 5a + 3/3-10a 2 -3/3 2 
-15a/?7 + 9<k*7 - 3a/35 + 10a 2 /35 + 15a 2 /3 7 + 24a/3 2 7 + 3f3-y5 + 4fij 2 5 - 95a 2 -/ 
-7a 3 06 - 5a 3 /?7 - 17a 2 /3 2 j - 4a(3 3 -/ + 5p 2 -/ 2 a + 35a 3 7 + 33a 2 /3 - 95a 2 + 16a/3 2 
-27a 3 /3 - 28a 2 /3 2 - 7/3 2 7 - 5/3 2 7 2 + 95a 3 - 4a/3 3 + 8a 4 /3 + 15a 3 /3 2 + 4a 2 /3 3 + 2/3 3 7 
+/3 3 7 2 - 35a 4 + 10a 3 - 5a 4 + /3 3 + a 5 , 

and 

P 52 = 1 - 25a7 + 6a/?7 - 3a 2 /?7 - 4a/3 2 7 + 5a 2 7 - 4a + 6a 2 - 2/3 - 4a 3 + a 4 + /3 2 + 9a/3 
-3/97 - 5a + 5 7 + 25a 2 - 12a 2 /3 - 4a/3 2 + 5a 3 /3 + 4a 2 /3 2 + 2/3 2 7 + /3 2 7 2 - 5a 3 . 

The following claim is proved using the Resolve function of the Mathematica system in Ap- 
pendix A. 4. 

Claim 25. P51 < 0, P52 > for all all (a,/3) € P and (7,5) in £/ie interior of (20). 



From Claim 25 we have that P52V-D — P51 > and hence the sign of P5 = P51 + P52V D is the 
same as the sign of DP 2 2 — P 2 X which is the same as the sign of 

, , n2 4(/3-5)(a- 7 ) 3 ((l-a) 2 5 + /3 2 (2a-l- 7 )) 2 

D - (P51/P52) 2 = — — J; — — > 0. 

-^52 

Hence P5 > which implies (33). □ 
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Proof of Lemma 18. Prom a + j3 + A(A — 2) a (3 < 1 we have 

(A - l) 2 < (1 - a)(l - p)/{atp) =: R n . (35) 

Lemma 24 implies that the coefficient of (A — 1) in (28) is positive and the coefficient of (A — l) 2 
in (28) is negative. Hence, also using (35), we have 

det(M) > lj.Rn [(-R 1+ R 2+ R 3 )(-R 1+ R 6 + R 7 )-Rl] + 

+ [(— -Rl — Rs — R9)(—Ri — Ra — R5) — R 2 ] I 

> — ^ — R1R2R11 — R±RqRu + -R1-R4 + RiR$ + 

+-R5-R9 + R3R7R11 + i?ii?5 + -R1-R9 — R\R%R\x — R\R-jR\\ + 
+R2R7R11 + R^RqRh ) . 



We will show 

— R1R2R11 — RiRqRii ~\~ R\R& + R\Rg > 0, (36) 

and 

-R5-R9 + R^RjR\\ + -R1-R5 + R±Rg — R1R3R11 — R1R7R11 + R2R7R11 + RsRqRu > 0. (37) 
To prove (36), note that 

— R1R2R11 — RiRqRh + R1R4 + RiR$ = 
D(l-a-/3) fl 1 (l-a)(l-/3) (l-a)(l-/3) 



l-a-/3-e-T/V7 5 a/3(l - 2a + 7) a/3(l - 2/3 + 5) 
The positivity of (36) follows from the following claim. 
Claim 26. 

1 1 (l-a)(l-/?) 
7 <5 a/3(l-2a + 7 ) a/3(l - 2/3 + <5) 

/or every (a, /3) € 7^3 D 7^a, (7><5) ^ the interior of (20), and (7, 5) in (23). 
Proof of Claim 26. 

-27<5 + <5a/3 + 7 a/3 - 4<5a/3 7 - 2<5a/3 2 + <5 2 a/3 - 25a 2 /3 + 4<5a 2 /3 2 - 25 2 a 2 /3 - 2 7 a/3 2 
-2 7 a 2 /3 + 4 7 a 2 /3 2 + 7 2 a/3 - 2 7 2 a/3 2 + 4 7< 5/3 - 2 7 <5/3 2 + 7< 5 2 /3 + 4 7< 5a + 7<5 2 a + 7 2 <5/3 
-27<5a 2 + 7 2 (5a - j5 2 - ~f 2 5 > 0, 

where the last inequality is proved using Resolve function of Mathematica system, see Appendix A. 5. 

□ 
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To show (37), we first note that 

R5R9 + R3R7R11 + -R1-R5 + R1R9 — R1R3R11 — R1R7R11 + R2R7R11 + R^RqRu 
_ / 4(l-/3-7-e) + (l-g)(l-/3) 2(g- 7 -g) (1 - a)(l - /3) 2(a - 7 - e) 



(a -7) (l-2a + 7)a/3 /8 - 8 (1 - 2/3 + S)a/3 a -7 

2(a-7-e) 2 /l 1\ {{I - a)(l - /3) 1 - /3 - 7 - 



+ T 



> 



(a — "f)(p — S) \t] e) \ a (3 a — j — i 

2y/D{a-i-i) ( (l-a)(l-/3) a - 7 - e / (1 - a)(l - /3) 1 - /3 - 7 - e 



/3 — 5 \ (1 — 2a + 7)a/3 (a — 7)77 \ a/3 a — 7 — e 

2-/D(a-7-e) / (1 - a)(l - /9) a - 7 - e / (1 - a)(l - /3) 1 - /3 - 7 - e 



a -7 \(1 - 2/3 + <5)a/3 (/3 - S)i \ a/3 a-j-i 

We will show 

(l-a)(l-/3) a -7-f^(i- a )(i-/3) 1 _ j g_ 7 _r, 



(1 — 2a + 7) a/3 (0 — 7)57 \ a/3 a — 7 — e 

for every (a, /3) E 7^3 D (7, 5) in the interior of (20), and (7, 5) in (23). Then by symmetry, we 

(l-a)(l-/3) a -7-^(l-a)(l_/3) i_0_ 7 _r, 



(1 - 2/3 + 5)a/3 08 -<f)e \ "Z 3 a-7-e 

(l-a)(l-/3) a-7-e /(l-a)(l-/3) 1 - /3 - 7 - e 



Let 

Py — ^ 

— 2a + 7) a/3 (0 — 7)57 \ a/3 a — 7 — e 

a/3?)(a — 7 )(1 — 2a + 7) 

= p 71 + p 72 Vd, 

where 

P 7 i = (3/2)a 2 + 4a/3 + (l/2)/3 2 -a<5 -7/3 + 7 5 + 1/2 -(3/2)a-/3 + ^ + (3/2)70/3 
+<5a/3 7 - 5a 2 /3 - (3/2) 7 a/3 2 - (l/2) 7 a 2 /3 - 7 5/3 - 7 5a - (l/2)a 3 - (9/2)a 2 /3 



and 



-(5/2)a/T + 7/3^ + a z 5 + (3/2)a J /3 + (5/2)cr/^ 



P 72 = (l/2) 7 a/3 - 1/2 + (l/2)/3 + a - (l/2)a 2 - (l/2)a/3 - (l/2)a 2 /3. 



The following claim is proved using the Resolve function of the Mathematica system in Ap- 
pendix A. 6. 

Claim 27. P71 > and P72 < for every (a,/3) 6 7^3 D R-a, (l,8) the interior of (20), and 
(7, 5) in (23). 

From Claim 27 we have that the following expression has the same sign as P 7 : 

(P /P \ 2 n _ til r T2 u 



pi 



where 



P 2 X - P 2 2 /J > 0. 



The last inequality is true for every every (a, /3) S R.3 D TZa, (7, 8) in the interior of (20), and (7, 5) 
in (23); it is proved using the Resolve function of the Mathematica system, see Appendix A. 7. □ 
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6.4 Proof of Lemma 19 

We first prove the following upper bound on (A — l) 2 . 

Proposition 28. For every A > 3, (a,f3) £ (7,(5) in the interior of (20), and (7,5) in (24), 
we have 

(A - l) 2 < (i-^X 1 -/ 3 ) < 1-P-l-e 
~ aj3 a — 7 — e 

Proof of Proposition 28. The first inequality is straightforward from a + P + A(A — 2) a/3 < 1. We 
next prove that (1 — a)(l — (3){a — 7 — e) < a/3(l — f3 — 7 — e). We have 

a/3(l - /3-7-e)-(l- a)(l - P)(a - 7 - e) 
= -(1- a- /3)(l-a- /3 + 2a/3 - V15) > 0, 



where the last inequality follows from 

D = (1 — a — P) 2 + 4(a — 7 )(/3 — <5) < (I — a — P + 2a/3) 2 . 

The last inequality is true when 7 > a 2 and 5 > p 2 (to see this note (1 — a — /3+2a/3) 2 — (1 — a — P) 2 = 
4{a - a 2 ){p - p 2 ) > 4(a - 7 )(/3 - <5)). □ 

Let i?io = (1 — /3 — 7 — e)/(a — 7 — e). Note that R7R10 = Rg, R3R10 = -R5. 
Lemma 24 implies that the coefficient of (A — 1) in (28) is positive and the coefficient of (A — l) 2 
in (28) is negative. Hence, also using Proposition 28, we have 

det(M) > 1(R 1Q ((-R 1+ R 2+ R 3 )(-R 1+ R 6 + R 7 )-Rl) 

+ (( — Ri — R% — -R9X— -Ri — R4 — R5) — R 2 )) j 
> — (-R 1 R 2 Rio - RiReRio + RiRa + RiRs + R5R9 + R3R7R10 + R2R9 + RM ■ 

We will prove that for every (a, (3) £ IZ3 D TZa., a, (3 < 1/2, (7, 5) in the interior of (20), and (7,(5) 
in (24), 

-R1-R4 — R1R2R10 + — {R5R9 + -R3-R7-R10) + R2R9 > 0, (38) 
and then by symmetry, we also have 

R\Rs — RiRqRio + - (RuRg + R3R7R10) + R^Rq > 0. 

The left-hand side of (38) can be written as the sum of the following two terms: 

(l-a-P)VD (\ 1 



1 — a — ft — e — 17 \ 7 1 — 2a + 7 
and 

v^D /2(1 -q)(1 -a-p-i) (1 - a - P) 



(1 — 2a + 7)57 \ a — 7 

We have (39) is non-negative, since 1 — 2a > 0. 



(39) 



(40) 
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For (40), we only need to prove 

l-a-P-e (l-a-/3) 2 
a — 7 i 

since 1 — 2a > 0. By (19), it is equivalent to prove 

l-a-/3-e-r)-(l-a-/3) 2 >0. (41) 

To prove (41), we have 

l-a-^-e-7?-(l-a-/3) 2 = -(l-a-/3) 2 -a-/3 + 7 + <5 + V~D, 

and 

D- (-(1 -a-p) 2 -a-/3 + 7 + 5) 2 = 2 7 + 25 - 2a 2 - 2/3 2 - 6a5 - 6 7 /3 + 2 7 5 - 2/35 - 5 2 

+6a/3 2 + 6a 2 /3 - 6a 2 /3 2 - 2a7 - 4a 3 /3 + 2a 2 7 + 2a 2 <5 
-4a/3 3 + 2/3 2 7 + 2/3 2 5 - 7 2 + 4a/35 + 4a/3 7 + 2a 3 
+2/3 3 - a 4 - /3 4 . 

The following claim is proved using the Resolve function of the Mathematica system in Ap- 
pendix A. 8. 

Claim 29. 

D - (-(1 - a - /3) 2 - a - /3 + 7 + S) 2 > 
/or every (a, 0) £lZ, a < 1/2, /3 < 1/2, (7, 5) m i/ie interior of (20), and (7,0") in (24)- 

since -(1 - a - p) 2 - a - p + 7 + 5 < 0, we can conclude that (41) is true. Hence, we complete 
the proof of Lemma 19. 

6.5 Proof of Lemma 20 

We have 

where the inequality follows from (20). 
We have 

(l_ a _ /3 )( a + /3 _ 7 _j) + ( / 3_ ( 5_ a + 7 )/D 

i?5 " i?7 " (a -7)0-*) • 

If/3 — 5 — a + 7>0, then i?5 > i?7. We may assume that /3 — 5 — a + 7<0. To show > R7, it 
is equivalent to prove that 



(1 - a - P)(a + p - 7 - 8) > {a - 7 - /3 + <5)VD. 

We have 

(l_ a _ j g)2 (a + j g_ 7 _g)2 4(a-7)(/3-g)(l-2/3-7 + g)(l-2a-5 + 7) 

(a -7-/3 + 5) 2 (a- 7 -/3 + <5) 2 

where the inequality follows from (20). 
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To show that R$ > AR% and Rg > ARj, it is sufficient to prove that 1 — (3 — 7 — i > 4(a — 7 — e), 
which follows from Proposition 28. 

We next prove U\ > 0. We first prove that Rg > 2Rq. We have 

where the first inequality follows from (20) and the last inequality follows from the fact that /3 < 
(3 - V5)/4. Hence ^ = R s - 2R 6 + R 9 - 2R 7 > 0. 
We next show: 

^ + R 9> 9R 2 . 

by the fact that i := |(1 + a — /3 — 27 — and (1 — a — /3) 2 = a/3, the goal is then to show, 

+ JL - £ ; > 



25 (3-5 l-2a + 7 

i.e., 

V^p 9(13-5) _ 3/3 



v / a/3 + 4(a- 7 )(/3-5) 1 - 2a + 7 25 

— a ' 

(i 



Since the above inequality is monotone with respect to 7, so let 7 = 7 m i n = - l^li 2 ^ , and the goal 



is to show: 

> 9(/3 ~ 8) _ M + 1/2 
Va^ + 4(a-7 min )(/3-<5) 1 - 2a + 7min 25 

or equivalently, 

1 9(/3-5) 3/3 



v / l + 4(l- 7min /a)(l-V/5) 1 " 2a + 7 min 25 
Replacing the term 7 m i n in the left-hand side, we get, 

1 > QW-S) _ M + 1/2 

^l + 4 (l-fiL^l)(l-J//3) l-2a + 7min 25 
Since a > 1/2 and 5 > /3 2 , replace them on the LHS, then our goal can be reduced to 

1 > 9 ^ M + i/2. 



V3-25//3 1 - 2a + 7 min 25 
Let X = 1 - 2a + 7m in = (irMilrgi+^zW) and A = So we are going to show 



25X > VA(185(f3 -5)- 3(3X + 5X), 



i.e., 

j2 v2 



45 2 X 2 > A(185((3 -5)- 3(3X + 5X) 2 , 

when 185(13 — 5) — 3/3X + 5X > 0. The rest of the proof can be checked by using Mathematica since 
they are all polynomial constraints, see Appendix A. 9. 
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A Mathematica Code 



Clear ["Global'*"] 

SetSystemOptions ["InequalitySolvingOptions" -> "CAD" -> True]; 
SetSystemOptions ["InequalitySolvingOptions" -> " Quadrat icQE" -> False]; 
SetSystemOptions ["InequalitySolvingOptions" -> "LinearQE" -> False]; 

RR = alpha > && beta > && alpha + beta < 1; 

RR3 = alpha > && beta > && alpha + beta + 3*alpha*beta < 1; 

Inter22 = gamma > && gamma < alpha && delta > && delta < beta && 

1 - 2*beta + delta - gamma > && 1 - 2*alpha + gamma - delta > 0; 
125 = gamma > && gamma < alpha" 2 && delta > && 

delta < beta"2 && (1 - alpha) ~2*delta + beta~2* (2*alpha - 1 - gamma) < && 

(1 - beta) ~2*gamma + alpha~2* (2*beta - 1 - delta) < 0; 
u26 = gamma > alpha~2 && gamma < alpha && delta > beta~2 && 

delta < beta && (1 - alpha) ~2*delta + beta~2* (2*alpha - 1 - gamma) > && 

(1 - beta) ~2*gamma + alpha~2* (2*beta - 1 - delta) > 0; 

DD = (1 - alpha - beta) "2 + 4* (alpha - gamma)* (beta - delta); 

epsilon = 1/2*(1 + alpha - beta - 2*gamma - SQ) ; 
eta = 1/2* (1 + beta - alpha - 2*delta - SQ) ; 

Wll = (alpha - gamma - epsilon)* epsilon* (1 - 2*alpha + gamma) / (eta* (alpha - gamma) 
W12 = (alpha - gamma)~2/((l - 2*alpha + gamma) *gamma) ; 

W21 = (beta - delta - eta)*eta*(l - 2*beta + delta) / (epsilon* (beta - delta)"2); 
W22 = (beta - delta) "2/ ((1 - 2*beta + delta) *delta) ; 

Rl = (1 - alpha - beta)/(l - alpha - beta - epsilon - eta); 
R2 = SQ/(1 - 2*alpha + gamma); 

R3 = 2* (alpha - gamma - epsilon) / (alpha - gamma); 
R4 = SQ/gamma; 

R5 = 2*(1 - beta - gamma - epsilon) /(alpha - gamma); 
R6 = SQ/(1 - 2*beta + delta); 

R7 = 2* (alpha - gamma - epsilon) /(beta - delta); 
R8 = SQ/delta; 

R9 = 2*(1 - beta - gamma - epsilon) /(beta - delta); 

RIO = (1 - beta - gamma - epsilon) / (alpha - gamma - epsilon); 

Rll = (1 - alpha)*(l - beta) / (alpha*beta) ; 
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A.l Mathematica Code for Proving Lemma 13 

W3 = Expand [Expand [Numerator [Wll] /2 - Denominator [Wll] /2] /. {SQ~2 -> DD>] ; 

Exponent [W3 , SQ] 

W31 = Coefficient [W3, SQ, 0] 

W32 = Coefficient [W3, SQ, 1] 

Resolve [Exists [{alpha, beta, gamma}, RR kk Inter22 kk W32 <= 0], Reals] 

W4 = Expand [Expand [Numerator [W21] /2 - Denominator [W21] /2] /. {SQ~2 -> DD>] ; 

Exponent [W4, SQ] 

W41 = Coefficients, SQ, 0] 

W42 = Coefficients, SQ, 1] 

Resolve [Exists [{alpha, beta, delta}, RR kk Inter22 kk W42 <= 0], Reals] 

W5 = FullSimplify[(W31/W32)~2 - DD] 
W6 = FullSimplify[(W41/W42)~2 - DD] 

A. 2 Mathematica Code for Proving Claim 22 

FullSimplify[((-Rl + R2 + R3) - ((beta - delta) * (-1/ (alpha - gamma - epsilon) 
- 1/eta + 1/epsilon) - SQ/epsilon + SQ/(1 - 2*alpha + gamma) + 
2*SQ/ (alpha - gamma))) /. {SQ -> Sqrt [DD] }] 

PI = Expand [FullSimplif y [ExpandAll [( (beta - delta) * (-1/ (alpha - gamma - epsilon) 
-1/eta + 1/epsilon) -SQ/epsilon + SQ/(1 - 2*alpha + gamma) + 
2*SQ/ (alpha - gamma) ) *epsilon*eta* (1 - 2*alpha + gamma) * (alpha - gamma)* 
(alpha - gamma - epsilon)] ]/. {SQ~2 -> DD, SQ"3 -> SQ*DD, SQ"4 -> DD~2}] ; 

Exponent [PI, SQ] 

Pll = Coefficient [PI, SQ, 0] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR kk Inter22 kk Pll <= 0], Reals] 
P12 = Coefficient [PI, SQ, 1] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR kk Inter22 && P12 >= 0], Reals] 

A. 3 Mathematica Code for Proving Claim 23 

P2 = FullSimplify[(Pll/P12)~2 - DD] ; 

P21 = -FullSimplif y [Numerator [P2]/( (beta - delta) ~2* (alpha - gamma) "3)] 
Resolve [Exists [{alpha, beta, gamma, delta}, RR && Inter22 && P21 <= 0], Reals] 

A. 4 Mathematica Code for Proving Claim 25 

P5=Expand [FullSimplif y [ 

Expand [((-Rl + R2 + R3)*(-R1 - R8 - R9) - Rl~2) *delta* (alpha - gamma)* 

epsilon* eta*(l - 2*alpha + gamma) * (alpha - gamma - epsilon)/ SQ]/. 

{SQ -> Sqrt [DD] }] ] / . {Sqrt [alpha" 2 + (-1 + beta) "2 + 

alpha (-2 + 6 beta - 4 delta) + 4 (-beta + delta) gamma] ->SQ} 

Exponent [P5, SQ] 

P51 = Coefficient [P5, SQ, 0] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR kk Inter22 && P51 >= 0]] 
P52 = Coefficient [P5, SQ, 1] 
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Resolve [Exists [{alpha, beta, gamma, delta}, RR kk Inter22 kk P52 <= 0]] 
FullSimplif y [DD - (P51/P52)~2] 

A. 5 Mathematica Code for Proving Claim 26 

Resolve [Exists [{alpha, beta, gamma, delta}, 
FullSimplif y [ 

Expand [(-R2*R11 - R6*R11 + R4 + R8) /SQ*alpha*beta*gamma* 

delta* (1 - 2*alpha + gamma) *(1 - 2*beta + delta)]] <= 
kk RR3 kk 125] , Reals] 

A. 6 Mathematica Code for Proving Claim 27 

P7 = FullSimplify[((l - alpha) *(1 - beta)/((l - 2*alpha + gamma) *alpha*beta) - 
(alpha - gamma - epsilon) /( (alpha - gamma) *eta) *( (1 - alpha)* 
(1 - beta)/(alpha*beta) - (1 - beta -gamma - epsilon) 
/(alpha - gamma - epsilon) )) *alpha*beta*eta* (alpha - gamma)* 
(1 - 2*alpha + gamma)] 

Exponent [P7, SQ] 

P71 = Coefficient [P7, SQ, 0] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR3 kk Inter22 kk P71 <= 0], Reals] 
P72 = Coefficient [P7, SQ, 1] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR3 kk Inter22 kk P72 >= 0], Reals] 

A. 7 Mathematica Code for Proving Lemma 18 
TP1 = Expand [P71~2 - P72~2*DD] 

Resolve [Exists [{alpha, beta, gamma, delta}, RR3 kk 125 kk TP1 <= 0], Reals] 
A. 8 Mathematica Code for Proving Claim 29 

TP2 = Expand [DD - (-(1 - alpha - beta)~2 - alpha - beta + gamma + delta) ~2] 
Resolve [Exists [{alpha, beta, gamma, delta}, RR kk alpha <= 1/2 kk beta <= 1/2 
kk u26 kk TP2 <= 0] , Reals] 

A. 9 Mathematica Code for Proving Lemma 20 

X = ((1 - 2*alpha)*(l - beta) "2 + alpha"2*(l - 2*beta + delta) )/(l - beta) ~2; 

A = (2*beta - 2*delta)/beta; 

Resolve [Exists [{alpha, beta, gamma, delta}, 

(1 - alpha - beta) "2 == alpha*beta kk alpha > 1/2 kk 
alpha < 1 kk beta > kk beta < 1 kk Inter22 kk u26 kk 
4*delta"2*X"2 <= A* (18*delta* (beta - delta) - 3*beta*X + delta*X)"2 kk 
18*delta*(beta - delta) - 3*beta*X + delta*X > 0], Reals] 
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